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Abstract 

We review recent progress on the phenomena of color superconductivity in high 
density quark matter. We begin with a brief overview of the unique aspects of physics 
near a Fermi surface and the imphcations for renormahzation group (RG) techniques. 
We next discuss the quahtative differences between asymptotic densities, where the 
effective couphng constant can be made arbitrarily small, and intermediate densities 
where quark matter is still strongly coupled. It is in the latter regime where RG 
techniques are particularly useful, in that they yield a generic description of possible 
behaviors without relying on an expansion in the strong coupling constant. Finally, 
we discuss aspects of the QCD groundstate at asymptotic densities, which can be 
determined in a systematic weak coupling expansion. 
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1 Introduction 



In this article I review recent progress on color superconductivity in high density QCD. It 
is remarkable how much more we know about QCD at high density compared to just two 
years ago, in 1998. The ultimate fate of baryonic matter at high density is a fundamental 
property of QCD, with implications for the astrophysics of neutron stars, as well as heavy 
ion collisions. More generally, we would like someday to understand the QCD phase diagram 
(figure (0)) in all its complexity as a function of temperature, chemical potential and quark 
masses. 

While lattice studies have been fundamental in determining the behavior of QCD at high 
temperature, technical difficulties arise in the application of Monte Carlo techniques once a 
chemical potential is introduced (essentially, the measure of integration is no longer positive 
definite, due to complex eigenvalues of the Dirac operator). The progress I describe does not 
rely on brute force techniques, but rather on physical insights associated with the presence 
of a Fermi surface. The special properties of physics near a Fermi surface are discussed in 
section 2. Remarkably, rigorous statements can be made in the limit of infinite density. 
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Figure 1: The QCD Phase Diagram as a function of temperature and density. 

There are two qualitatively different regimes to be addressed. The regime most likely to 
be realized in nature (i.e. in the cores of neutron stars or in heavy ion collisions) is that 
of intermediate density, where the average distance between quarks is still of order a Fermi 
and the gluons exchanged are rather soft. This regime is strongly coupled, and while we 
can make qualitative statements about the nature of the ground state using renormalization 
group techniques, we cannot perform quantitative calculations. The intermediate density 
case is discussed in section 2. 
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At extremely high densities the effective couphng is small, leading to a weakly coupled 
liquid of quarks. The only low energy excitations in this liquid are quasi-particles and quasi- 
holes representing fluctuations near the Fermi surface. There are rather general arguments 
(see below) which suggest that any attractive interaction, no matter how small, can lead to 
Cooper pairing in the presence of a Fermi surface. In the case of QCD, this attractive inter- 
action is provided by gluon exchange in the 3 channel. Thus, we expect to find condensation 
of Cooper pairs of quarks in the high density limit. 

Actually, even in the high density limit, there are some subtleties which make the problem 
less than straightforward. Most importantly, although the effective coupling is weak, the 
exchange of magnetic gluons leads to a long range interaction due to the absence of a magnetic 
screening mass. Quarks carry electric, rather than magnetic, color charge, and hence are 
better at screening the timelike component (Aq) of the gluon field than the spacelike (Ai). 
It is only non-perturbative effects which can screen color- magnetic fluctuations, and these 
are nearly absent at weak coupling. An understanding of dynamic screening due to Landau 
damping is necessary to control the high density calculations. These issues are discussed in 
section 4. 

The order parameter for color superconductivity is 



where C is the charge conjugation operator and F is a matrix in color, flavor and Dirac space. 
As we will see below, determining the precise form of F requires some work; it is particularly 
complicated in the case of three flavors. However, as we will discuss in section 5, in the 
weak coupling limit the true groundstate can be determined in a controlled approximation. 
Here we will simply note that single gluon exchange between two quark quasiparticles can 
be decomposed into an attractive 3 channel and a repulsive 6 channel. Thus, at the most 
naive level we expect an anti-triplet condensate, which breaks SU(3) —>■ SU(2). The main 
effect of this condensate is that it leads to the Higgs phenomena for at least some subset of 
the gluons, or equivalently to the Meissner effect and screening of some subset of the color 
magnetic fields. The phenomenological implications of color superconductivity are still not 
well understood and will be the subject of investigation for some time to come. 

Let me close this introduction with a historical note. The idea that quark matter might 
be a color superconductor is quite old The original insight was based on the existence 

of the attractive 3 channel and an analogy with ordinary superconductors. Recent interest in 
the problem was rekindled by the work of two groups that considered diquark condensation 
due to instanton- mediated interactions [|], predicting gaps as large as ~ 100 MeV. These 
calculations, while uncontrolled, are quite suggestive, and led to the recent progress on the 
subject. It is often claimed that early investigations predicted tiny gaps, at most of order a 
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few MeV. While this may have been the consensus among the few theorists who had actually 
worked on the problem, it is actually an unfair characterization of Bailin and Love's results 
0. A value of the strong coupling large enough to justify the instanton liquid picture of [|] 
also yields a large gap when substituted in Bailin and Love's results. After all, instantons 
are suppressed by an exponential factor exp(— 27r/Q;s). Bailin and Love merely suffered from 
the good taste not to extrapolate their results to large values of a^! 



2 Physics Near a Fermi Surface and the Renormahza- 
tion Group 

Important simplifications arise in the study of cold, dense matter due to the existence of a 
Fermi surface, which in relativistic systems is likely to be rotationally invariant. The energy 
of a low-energy excitation (a quasi-particle or -hole) is then independent of the orientation 
of its momentum p, and only depends onp — fi, where p = |p| and /i is the chemical potential 
or Fermi energy. (Here, for simplicity, we will always work with massless quarks.) This 
leads to a kind of dimensional reduction, so that physics near a Fermi surface is effectively 
1+1 dimensional. In particular, arbitrarily weak interactions can lead to non-perturbative 
phenomena like pair formation. 

The renormalization group approach is particularly useful here - we integrate out the 
modes far from the Fermi surface, leaving only the low-energy quasi-particle and -hole states 
that are involved in the interesting physics. These excitations might in principle be related 
to the original quarks in a complicated way, but on quite general grounds must be described 
by an effective action of the form 

Seff = J dt d?p {tdt - (e(p) -eF))ij + S^nu (2) 

where Sint contains higher dimensional, local quasi-particle operators. Strictly speaking, 
this form of the effective action is only valid for models in which the original interactions 
were local (short ranged). While appropriate for QCD at intermediate densities 0, where 
non-perturbative effects are expected to generate screening of magnetic gluons, it must be 



modified at weak coupling where magnetic fluctuations are long ranged |jTO|, [TT[]. However, 
it is the only technique I know of from which we can obtain robust information about the 
strongly coupled region of the phase diagram. Below I review the results of this analysis, 
and defer a discussion of the weak coupling phase until the following section. 

It can be shown |^, |^ using simple classical scaling arguments that all interactions are 
irrelevant except for the Cooper pairing interaction (scattering of quasi-particles at opposide 
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sides of the Fermi surface: pi ~ —P2) and stricly colinear scattering: pi — P2, which can 
lead to the Overhauser effect (chiral waves) at large- A*",, |p. Both of these interactions are 
classically marginal, so quantum corrections determine their evolution. Here we restrict 
ourselves to local Cooper pairing operators which are invariant under the full SU{3)l x 
SU{3)r X U{1)a chiral symmetry: 

OIl = (^l7o^l)', OIj, = {^PLlo^L)i^pRlo^R) (3) 

These come in both color symmetric (3) and antisymmetric (sextet) combinations. More 
general operators with different flavor or Dirac structures can be reduced to linear com- 
binations of the basic ones (§), using parity and Fierz rearrangements. This analysis can 
be extended to operators (such as those induced by instantons) that break the anomalous 
U{1)a symmetry p, 0, yielding a very robust characterization of QCD even at intermediate 
densities and strong coupling. We will not discuss the details of this more general analysis 
here, but the results (given reasonable assumptions about the signs and magnitudes of the 
interactions) are qualitatively similar 

The RG evolution of the operators in (|]) is determined by quark-quark scattering near 
the Fermi surface. A bubble graph with four-quark vertices Fi and F2 and external quark 
lines satisfying the Cooper pairing kinematics yields 



G1G2I {T^),^{T 



1 k'k 



ilo)ij{lo)kl + ^{l)ij{l)kl 



r2)jA^2)u' (4) 



Here / = g;^/x^ log(A//j/Af/y), where [AfR,Ai/v] are the upper and lower cutoffs of the 
momentum shell integrated out. We define the density of states on the Fermi surface to be 
A^ = /i^/(27r^) (in weak coupling) and t = log(A//j/A[/y). The RG flow does not mix LL 
and LR operators, nor different color channels. We obtain the following RG equations 

— It — - ^ ' 



dt 

d{G^^ + 3Gf ^) 



dt 

d{GQ^ - Gf^) _ 2N ^^^^ r^LR\2 
dt 



NiGi^^ - 3Gf^)^ (6) 
0, (7) 
(Go^^ - Gf 



The linear combination G^ = Gq^ + Gf^ reaches its Landau pole first, governed by the 
equation 

l + (Ar/3)a(0)t ■ 
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In general, interactions which are attractive at the matching scale, G'*(0) > 0, will grow 
during the evolution, and reach a Landau pole at the scale t^, = — 3/(A^G'*(0)). The corre- 
sponding energy scale is 

A,„ = A„vexp(--^). (10) 

which agrees with the usual BCS result. Repulsive interactions, corresponding to negative 
initial values of the coupling, become weaker near the Fermi surface. 




Figure 2: Excitations on opposite sides of the Fermi surface. 



3 High Density Limit 

In the high density limit the typical momentum transfer between quarks is large, and there- 
fore the effective coupling is small. The properties of this phase can be deduced in a sys- 
tematic, weak coupling expansion. However, there is one technical problem that must be 
solved having to do with soft magnetic gluons. In the renormalization of the Cooper pairing 
interaction there is a region of phase space where the incoming quarks are only slightly de- 
flected by the gluon exchange. This leads to an IR divergence unless arbitrarily soft gluons 
are screened in some way. While gluons acquire a perturbative electric mass at high density, 
it can be shown that to all orders in perturbation theory no magnetic mass is generated 0. 
The only hope is that Landau damping - a form of dynamic screening affecting the spacelike 
gluons - is enough to control this IR problem. That this is so was first pointed out by Son 



lOl , who went on to deduce the foUowing behavior for the diquark gap at weak coupling 

V2g 



A~;.,-^exp(-^) . (11) 
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This result has since been confirmed by RG methods [|TTI] as well as Schwinger-Dyson tech- 



niques [12]. In this section I discuss some of the details of these calculations, concentrating 



on the the RG approach. 

The magnetic gluon propagator, including vacuum polarization effects from virtual quarks, 
has the form (go « 

DUqo,q) = I I • (12) 

Strictly speaking D^^ is gauge dependent, but in our leading order calculations the propaga- 
tor always appears contracted with gamma matrices next to nearly on-shell external quark 
lines. Thus the gauge dependent parts are higher order in g due to the equations of motion. 
The effect of Landau damping is to cut off the small-q divergence in (|12]) at g ~ g^^m^^, 

2 2 

where = -^/^^ is the Debye screening mass. A common feature of both the Schwinger- 
Dyson and RG calculations in the weak coupling region is loop integrals dominated by energy 
transfers of order go ~ A, and hence momentum transfers of order g* ~ g* can be 

made as large as desired by going to high density. 

The main technical problem in the RG approach is long range magnetic interactions, or 
equivalently the presence of soft gluons. At no point can the theory be completely described 
by quarks with purely local interactions as in (^. The effective Lagrangian contains both 
quark and gluon excitations (with energies below the cutoff A) and local interactions resulting 
from integration of higher energy shells. This modifies the form of the RG equations obtained 



11|, so that equations (O) and (BI) become 



'^^"'^t^^-'^ = -f (G- + G-)^ - 1^ , (13) 
= -N{Gl;^ - 3Gf-)^ + ^ . (14) 



The solution of these RG equations leads to a Landau pole in the dominant 3, LL and 
RR channels given by ([111) . It is worth commenting on the angular momentum of the 
condensate. The RG equations can be derived for general values of angular momentum 1. 
Naively interpreted, the results suggest that condensates might occur in higher 1 channels, 
leading to the breaking of 0(3) rotational invariance. A more detailed gap equation analysis 



111 shows that this is not the large s-wave gap suppresses the formation of p-wave 

and higher 1 gaps. The literature is somewhat confused on this important issue. The papers 
in |jl3| address the issue of rotational invariance and do not agree on the size of higher 1 gaps. 
However, neither paper addresses the interplay of s-wave and higher 1 gaps, which is studied 
in O. 
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4 The QCD Groundstate at High Density 



In this section I describe the vacuum energy analysis necessary to determine the groundstate 
of QCD at high density |jT^. Neither the RG nor Schwinger-Dyson analyses are sufficient to 



specify the actual groundstate. Strictly speaking, the former only reveals the energy scale 
and quantum numbers of the pairing instability, while the latter only identifies extrema 
of the vacuum energy. As we shall see, there are additional subtleties which can only be 
resolved by consideration of energetics. 

First, let us consider the case of 2 massless flavors. Because the condensate occurs 
between pairs of either left (LL) or right (RR) handed quarks in the J=L=S=0 channel ||11|| , 
and the 3 color channel is antisymmetric, the quarks must pair in the isospin singlet (ud 
- du) flavor channel. However, even in this case there is a subtlety, as the relative color 
orientations of the LL and RR condensates are not determined by the usual leading order 
analysis. A misalignment of these condensates violates parity, and further breaks the gauge 
group beyond SU(3)c SU(2)c. An analysis of the Meissner effect is necessary to determine 
the relative orientation [|l^, and the effect is higher order in g. There are thus a number of 
unstable configurations of only slightly higher energy with different color-flavor orientations 
(and hence different symmetry breaking patterns), leading to the possibility of disorienting 
the diquark condensate (see figure (0)). 





L-L R-R 

Figure 3: Color disorientation of LL and RR condensates. 





ud-du us-su ds-sd 

Figure 4: CFL condensation: each flavor channel orients chooses an orthogonal color orien- 
tation. 

The generalization to three flavors is far from straightforward. Again, one can show that 
the condensate must occur in the J=L=S=0 and color 3 channel. (The sextet condensate is 



suppressed in the weak coupling limit [0] and I do not discuss it here.) The Pauli principle 
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then requires that the flavor structure again be antisymmetric ~ (g^gj — qjqi), for quarks of 
flavor Thus, one can have combinations of condensates which are in the 3 of both color 
and flavor SU(3)l or SU(3)r. Due to the chirality preserving nature of perturbative gluon 
exchange, there is no mixing of LL and RR condensates, which form independently. One can 
immediately see that there are a number of possibilities. For example, the condensates for 
the three flavors and both chiralities might all align in color space, leading to an SU(3)c — > 
SU(2)c breaking pattern. A more complicated condensate has been proposed [|l^, |l^ called 
Color Flavor Locking (CFL), in which the 3 color orientations are "locked" to the 3 flavor 
orientation. In figure (H) we give a simple picture of CFL condensation. 

To determine the nature of the energy surface governing the various color-flavor orienta- 
tions of the condensate, we can begin by characterizing the color-flavor configuration space 
of condensates. We consider the ansatz 

Aab L,R AC L,R Mbc /i r\ 

ij — ^ ^ijk , U'^j 

where a,b are color and i,j flavor indices. L and R denote pairing between pairs of left and 
right handed quarks, respectively. Under color and flavor A transforms as 

U^A^V^ , (16) 

where Uc is an element of SU(3)c and of SU(3)l. A similar equation holds for A^. It is 
always possible to diagonalize A^ by appropriate choice of Uc and V^: 



A' 



/ a 0\ 
6 

V c y 



(17) 



Generically, there does not exist a which diagonalizes A^ in this basis. In the CFL 
case, where the diagonalized A^ is proportional to the identity, a = 6 = c, it is easy to 
show that one can choose such that A^ = ±A^. These two configurations are related 
by a U{1)a rotation (see section 3). Hence, they are degenerate in the high density limit 
where gluon exchange dominates. Instanton effects, important at intermediate density, favor 
A^ = A^. Note that parity, if unbroken, requires A^ = A^, and hence implies simultaneous 
diagonalizability. 

In 0] we considered the potential vacua parametrized by a,b,c. First, we use the Dyson- 
Schwinger (gap) equation to determine which of these configurations are energy extrema. 
Next, we computed the energies of the extrema to determine the true groundstate. A similar 



analysis has been carried out by Schafer and Wilczek |]l6l in the approximation where gluon 



interactions are replaced by local four fermion interactions. They concluded that the CFL 
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vacuum had the lowest energy. In our analysis, which I summarize below, we included 
the gluons in the analysis, introducing long range color-magnetic fluctuations (controlled 
by Landau damping) and Meissner screening into the gap equation and vacuum energy 
calculations. 

At asymptoticaly high densities (weak coupling) the diagrams (a)-(c) in figure (^ give 
the leading approximation to the effective action. Note that in these diagrams the quark 
propagators include the diquark condensate (see (^0|) below), and the gluon propagators 
include Landau damping, but not the Meissner effect. The latter arises from the condensate- 
dependence of quark loops in diagrams (c) and (d). The resulting gap equation (figure (|)), 
with condensate shown explicitly at lowest order in A) is given by 



D'^^ is the gluon propagator, including the effects of Landau damping and Debye screening 
(we assume Feynman gauge throughout). 

We will restrict the color group structure in the gap equation to the attractive anti- 
symmetric 3 channel, which projects out the anti-symmetric part of S{k) in color space in 
the gap equation. Here S is the fermion propagator for the spinor with i a flavor 

index and a a color index. 




(18) 



where 




(19) 




Figure 5: Vacuum energy diagrams. 





Figure 6: Schwinger- Dyson equations. 
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For the three flavor case S can be written exphcitly as an 18 x 18 matrix in color flavor 
space. The inverse propagator may be written 

\ + 7oAt7o 
A 



(20) 



where = fijo- A is a 9 x 9 matrix which for the ansatz (|T 



A 






c 










— c 















-b 








— c 










c 







a 











—a 







-b 















—a 






takes the form 

b \ 



a 







(21) 



Because we are deahng with a diquark condensate the non-trivial part of the gap equation 
involves the lower left 9x9 block. We will refer to this sub-block of the propagator S as S2i- 

For a particular ansatz A to be a solution to the gap equation we require that the 
color antisymmetric part of S2i{k) (corresponding to the 3 channel) be proportional 
in color-flavor space to the off-diagonal 21 submatrix of S~^{q) — SQ^{q), or A(g), which 



appears on the LHS of the gap equation (see fl^ for more discussion of this point). 



The propagator may be found by inverting the sparse matrix in (|20| ) using Mathematica. 
Only three ansatze satisfy our condition: a = b = c; a = b, c = 0; b = c = 0. We refer to these 
solutions as (111) (color-flavor locking), (110) (3^0 breaking) and (100) (3^2 breaking) 
respectively. For these ansatze the color antisymmetric part of T^S2i{k)T^ has the form of 



a constant multiplying the matrix form (^) with a, 6, c set to or 1 as is appropriate for 
the ansatz. After contour integration over 1, we flnd the following gap kernels 

2 A 1 A 

3' 



(111) 
(110) 
(100) 



A 



A2 



A 



A2 



+ 



ki^ 
A 



4A2 



2jki 



2A2 



A2 



(22) 



These kernels are to be substituted in the following gap equation, which we obtain under 
the approximation go << (We also neglected the anti-particle contributions, which are 
suppressed by powers of 1/yU.) 
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2 . . / 3 



— ^ COS 9 



A(po) = I dqo I dcosO { , , ' ,,,,,, (23) 



127r2y ^"J \i-cose + {G + ipo-qo)^)/i2i2 



1 

2 ' 2 



i COS 6* 



l-cose+{F+{po-qoy)/{2l2^ 



Kiqo), 



where F and G represent the medium effects on the electric and magnetic gluons, and K{qQ) 
is one of the gap kernels from (|2^) . The Meissner effect makes an additional contribution 
to G beyond that of Landau damping. In \TZ\ we evaluated the gluon vacuum polarization 



P^iij{qo,q) in the presence of a diquark condensate. (A more detailed computation of the 
Meissner effect is given by Rischke [jl^, with similar results.) The additional Meissner 



screening is given by 6G = \VjjPij , where Vj^ = (dij — qicjj^ is the transverse projection 
operator. At low momenta, qo, g ~ A, 6G{qo, q) is of order the Debye mass mo ^ gfJ', while 
at larger energy or momenta the effect is suppressed by a power of ^ or ^. We limited 
ourselves to an estimate of the size of the Meissner effect on the gap solutions. To this end, 
we used the following approximation for 6G: 

5G{q,,q)^ml ^° , (24) 

+ ql + 

where Aq is the maximum value of the function A(A;o, k). Note we did not introduce any 
color structure in 5G\ all gluons experience the same magnetic screening. While this is a 
crude approximation, it gives the rough size of the Meissner effect on A. 

We solved the gap equations for all three gap kernels using this form of the Meissner 
effect, and the results are shown in figure (0) for the case of /i = 400 MeV. The effect is to 
decrease the size of the condensate but it is a small perturbation on the solutions obtained 
without the Meissner effect. 

To determine which of the above gaps is the true minimum energy state we must cal- 
culate the vacuum energy, which receives contributions from vacuum to vacuum loops of 
both quarks and gluons (figure 1). We start with the CJT effective potential |0, which 
upon extremization wrt appropriate propagators and vertices leads to the Schwinger-Dyson 
equations. The fermion equation is the gap equation given above, while the gluon equation 
reproduces Landau damping. We wish to compare energies corresponding to our three solu- 
tions to determine which one is the true vacuum (the difference in energies V will be gauge 
invariant, whereas actual values are not). It is easy to show that the value of the effective 
potential evaluated on the gap solution is given by: 

V = -^J-0^,tT\nS{p)/So{p). (25) 
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Figure 7: Gap solutions for ^ = 400 MeV, with and without Meissner effect. 



Diagramatically, this is equivalent to the graph of figure 1(a) when evaluated on the gap 
solution. 

The fermion loops are most easily calculated by going to a basis where SqS'^ is diagonal 
in color-flavor space. Note that the gap matrix A has non-trivial Dirac structure that must 
be accounted for |T^: A = A175P+ + A275P-- , where P± are particle and anti-particle 
projectors. Our analysis has been restricted to the particle gap function Ai. The anti- 
particle gap function A2 has its support near fco ~ 2/i, and its contribution to the vacuum 
energy is suppressed. There are 18 eigenvalues, which occur in 9 pairs. The product of each 
pair is of the form 



- 1 + a- , , 

where a is an integer. For our three cases we obtain the following sets of eigenvalues: 



(26) 



(111) - 


8 X 


{a 


= 1}, 


1 X {a = 4} 


(110) - 


4 X 


{a 


= 1}, 


2 X {a = 2} 


(100) - 


4 X 


{a 


= 1} 





(27) 



The binding energy is of order 



- I d kdko In 



1 + a 



A\ko,k) 
kl + ik- 



(28) 
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~ —a /i^Aq 



(29) 



where Aq is the maximum value of the gap function A(/co, k), which has rather broad support 
in both energy and momentum space away from the Fermi surface. A more precise answer 
than ( ^8]) requires numerical evaluation, but it is clear that the result scales with a and has 
only a weak (logarithmic) dependence on the variations in the shape of A(A;o, k). Substituting 
our numerical results for the gaps in the three cases, it is easy to establish that 



We find that the CFL vacua remains the lowest energy state, at least at asymptotically 
high densities where the calculation is reliable. The Meissner effect is a small correction to 
the vacuum energy at asymptotic densities. Configurations which satisfy the gap equations 
but are not the global minimum of energy are presumably saddlepoints, since they are 
continuously connected to the CFL vacuum via color and flavor rotations. 

5 Conclusions 

In this contribution I have tried to summarize some important progress of the last two years 
on the theory of cold, dense quark matter. Due to space limitations, I was not able to discuss 
a number of important issues, such as the low energy effective Lagrangian, continuity of 
hadronic and quark phases and more phenomenological studies. I list some of the important 
papers in ||20[| . 

I thank my collaborators N. Evans, J. Hormuzdiar and M. Schwetz and my colleagues D. 
Hong, R. Pisarski, K. Rajagopal, M. Rho, D. Rischke and T. Schafer for contributing to my 
understanding of this subject. My research was supported under DOE contract DE-FG06- 
85ER40224 and by a JSPS visiting fellowship. 
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